The Plant Journal

Whether Photosystem I (PSI) core subunits accumulate prior to light exposure in developing angiosperm seedlings remains unresolved, with conflicting reports across species. Here, we investigated the presence and membrane colocalization of the PSI core subunit PsaA in etioplasts of dark-grown angiosperms representing dicot and monocot species. Immunoblotting showed that PsaA accumulates in etioplasts of all three dicot species examined (pea, Arabidopsis, and runner bean), whereas in the monocot oat it was detected only after prolonged etiolation, at substantially lower levels and with an anomalously high apparent molecular weight. Blue-native PAGE analysis reveals that a fraction of PsaA co-migrates with LPOR, PsaB, FNR, and chlorophyll synthase, suggesting co-localization within a shared membrane microdomain rather than stable complex formation. The thylakoid insertase Alb3 was more abundant in dicot etioplasts, consistent with a potential role in the early integration of PsaA into the membrane. Upon illumination, pea reached PSI functionality faster than oat, with P700 oxidation detectable 30 min earlier, linking the dark accumulation of PsaA to an accelerated photosynthetic onset. These findings demonstrate light-independent accumulation of a PSI core subunit in a species-dependent manner and point to early steps in PSI biogenesis that precede full photosynthetic complex assembly. Highlight Contrary to prevailing models, a Photosystem I core subunit PsaA accumulates in dark-grown angiosperm seedlings before light exposure, revealing light-independent early steps in photosynthetic complex biogenesis.

The genus Rhynchospora Vahl (beak-sedges) comprises approximately 381 accepted species with a worldwide distribution, all of which possess holocentric chromosomes, where centromeric activity is distributed almost along the entire chromosome. Despite the recent advances, the mechanisms governing the dynamics of meiotic recombination in holocentric plants remain poorly understood. Here, we developed haplotype-specific oligo-FISH probes for chromosomes 1, 2, and 3 based on a haplotype-phased genome assembly of Rhynchospora breviuscula (n = 5), enabling homolog-specific chromosome painting. Each probe set was labelled with a distinct fluorophore and hybridised in situ to metaphase chromosomes of the reference plant and seven F1 individuals derived from self-crossed reference plants. This approach allowed the unambiguous discrimination of homologous haplotypes and the indirect visualisation of crossover (CO) events in recombined chromosomes. We observed that recombination events were predominantly located in terminal chromosomal regions, consistent across individuals. These results corroborate previous findings from single-cell recombination mapping and provide independent cytological validation of the recombination landscape in this species. Our study establishes haplotype-specific chromosome painting as a robust tool for high-resolution mapping of meiotic recombination in holocentric plants across generations. Furthermore, these probes provided a foundation for future investigations into inverted meiosis, a mechanism characterized by an alternative pattern of chromosome segregation in holocentric species.

In Asteraceae species, the capitulum is a compact inflorescence, featuring a characteristic reproductive structure. Despite the identification of a few key regulatory factors, the transcriptome-level information on the developing capitulum remains limited. Here, we applied single-cell and spatial transcriptome sequencing to investigate the developing Gerbera hybridas capitulum during floret differentiation. We obtained a transcriptomics atlas encompassing different stages of the Gerbera capitulum and analyzed the cellular and spatial dynamics of gene expression. Using marker gene expression and GO enrichment of cluster-specific DEGs, we annotated putative cell types and described changes in gene expression across sampled stages, potentially associated with ongoing developmental processes. We detected activity of previously undescribed MADS-box genes and defined their spatial expression patterns. Notably, the MADS-box gene GAGL12 was found to be enriched in the putative capitulum phloem cells. The GAGL12 protein was shown in yeast two-hybrid assays to interact with several other MADS-domain proteins with hypothesized functions in vasculature development, and further detailed in silico analyses supported a candidate role in the development of capitulum vasculature. Altogether, we provide integrative and dynamic transcriptomic insight into capitulum and floret development and lay a basis for future functional studies of the control and development of this intriguing reproductive structure.

Auxin is a key phytohormone that regulates all aspects of plant growth, development, and environmental responses, making the precise analysis of its distribution and signaling essential for understanding plant adaptation and physiological processes. However, despite the agricultural importance of oilseed rape (Brassica napus), the lack of robust, species-specific molecular tools limits detailed studies of hormone signaling in this crop. Here, we developed and characterized reporter systems for the sensitive visualization and quantification of auxin distribution and signaling in B. napus. The DR5cc auxin signaling reporter and a novel synthetic auxin-responsive reporter, BIP3, assembled from promoter fragments of three oilseed rape IAA genes, were generated to drive GUS expression. In hairy roots, both reporters showed auxin-responsive expression in the root apical meristem that became broader after auxin treatment. In transgenic seedlings, flowers at anthesis, and 12-day-old embryos, DR5cc exhibited a more defined expression pattern than BIP3. To monitor real-time auxin dynamics under abiotic stress, DR5cc fluorescent reporters were employed in hairy roots. Mannitol and NaCl treatments induced a time-dependent increase in fluorescence, peaking at 6-12 h before returning to basal levels after 24 h. Furthermore, dual-reporter assays enabled simultaneous monitoring of auxin and cytokinin signaling, revealing distinct hormone-specific spatial responses in hairy roots. Finally, we established a quantitative DII (qDII) reporter system using degron domains from B. napus Aux/IAA proteins, providing a high-resolution quantitative readout of auxin depletion. Together, these reporter systems enable spatial, temporal, and quantitative analyses of auxin dynamics during development and stress adaptation in oilseed rape.

Leaf senescence directly affects lettuce quality and postharvest shelf life, but the regulatory roles of miR319-targeted and non-target CIN-TCP transcription factors remain unclear. Here, we examined whether the miR319-TCP module controls lettuce leaf senescence through separable genetic branches. MIR319 overexpression delayed dark-induced senescence, whereas STTM-mediated miR319 suppression accelerated chlorophyll loss, photosynthetic decline, and senescence-marker activation. Disruption of the miR319-targeted gene LsTCP4 phenocopied MIR319 overexpression, supporting LsTCP4 as a pro-senescence factor downstream of miR319. We further found that the miR319 non-target CIN gene LsTCP17 also promoted senescence, as tcp17 leaves retained more chlorophyll than wild type during dark treatment. Genetic combinations showed that tcp17 enhanced chlorophyll retention in the OX319 background and partially rescued the accelerated senescence phenotype of S319, indicating that LsTCP17 acts through a route separable from the miR319-targeted branch. Together, these results reveal a split CIN-TCP architecture in which miR319-targeted LsTCP4 and non-target LsTCP17 provide parallel pro-senescence inputs, offering a genetic framework for targeted improvement of lettuce quality.

O_LIHybridisation is frequently associated with plant invasions; however, its consequences for genome organisation and chromosome evolution remain poorly understood in invasive species. We investigated the extent of hybridisation in the invasive Carpobrotus edulis--acinaciformis hybrid complex and determined the cytogenomic contribution of parental species in hybrid accessions. C_LIO_LIWe combined whole-genome sequencing, population genomic analyses, genome size estimation, repeatome characterisation, chromosome counting and fluorescence in situ hybridisation to compare parental species and hybrid accessions from South Africa and the Mediterranean Basin. C_LIO_LIPopulation genomic analyses revealed widespread hybridisation and introgression, with most invasive accessions showing admixed ancestries. Pattersons D-statistic supported asymmetric allele sharing towards C. edulis. Hybrid accessions displayed genome sizes indistinguishable from C. edulis, whereas C. acinaciformis possessed significantly larger genomes. Repeatome analyses identified marked differences in repetitive DNA composition, particularly in satellite DNA abundance and chromosomal distribution. A newly identified satellite repeat (CarpoSat) showed contrasting chromosomal patterns between parental species, whereas hybrids resembled C. edulis satellite pattern. C_LIO_LIOur results demonstrate that Carpobrotus hybrid accessions are a swarm of later-generation hybrids and backcrosses showing a strong bias towards C. edulis, indicating asymmetric introgression. These findings highlight the value of integrating cytogenetic and genomic approaches to understand genome evolution in invasive hybrid complexes. C_LI

The receptor kinase FERONIA (FER) perceives endogenous RAPID ALKALINIZATION FACTOR (RALF) peptides and regulates a plethora of plant physiological processes, including immunity. RALF peptides also bind to LEUCINE-RICH REPEAT EXTENSIN (LRX) proteins as structural components of the cell wall. We recently showed that the FER-RALF-LRX pathway supports colonization by the obligate biotrophic fungal pathogen Erysiphe cruciferarum (Ecr), a member of the powdery mildew species complex that infects Arabidopsis. Genetic disruption of the pathway primarily affects conidiation of the fungus, raising the question of effects on fungal nutrition. To get further insight into the underlying mechanisms, we performed RNA sequencing (RNAseq) to identify differential transcriptional responses of FER-RALF-LRX pathway mutants upon Ecr infection. Surprisingly, our results revealed that pathway disruption has a limited impact on the overall transcriptional changes upon fungal infection. However, consistent with previous reports, FER-RALF-LRX pathway mutants show changes in basal expression of a plethora of genes, mainly associated with cell wall metabolism, jasmonic acid signalling, amino acid biosynthesis and secondary metabolism. Many of these genes are regulated by Ecr infection across genotypes, too. This raises the question whether these are relevant pathway components for powdery mildew host establishment downstream of the FER-RALF-LRX module. In summary, our data reveals new insights into FER-RALF-LRX-dependent responses that may support host susceptibility to biotrophic plant pathogens.

BackgroundStructural variants (SVs) are large insertions or deletions of DNA sequences. While less numerous than single nucleotide polymorphisms, SVs often account for a greater proportion of nucleotide differences between genomes. Their size and frequent association with repetitive sequences has historically hindered their detection, which has limited the ability to associate this variation with molecular and phenotypic trait variation. While some SVs have been linked to observable traits, it remains unclear whether such effects are rare or broadly distributed across the genome. ResultsTo test for genome-wide relationships between SVs and gene expression, we analyzed genome assemblies and transcriptomic data from 10 tissues across 26 diverse maize inbred lines. We identified SVs amongst these lines and examined variants located within the 1kb promoter region upstream of genes. Thousands of genes showed expression differences associated with promoter SVs, often in a tissue-specific manner. One common feature of these SVs was the presence of transposable element sequences. LTR retrotransposons were enriched amongst promoter SVs associated with differential expression and often reduced expression of the nearby gene. Despite widespread expression changes, we found no enrichment for specific biological functions or pathways among affected genes. ConclusionsOur findings indicate that extant TE-mediated promoter SVs play a significant role in shaping gene expression patterns across the maize genome. However, their phenotypic effects appear limited or context-dependent, suggesting that many variants may have minimal impact outside specific developmental stages or environmental conditions.

Elongating rice internodes present a developmental gradient from dividing meristem to mature cells, providing an elegant pseudo-time course for study of plant vegetative development. We tested the hypothesis that DNA methylation regulates gene expression during rice internode development by integrating RNA-seq and bisulfite DNA sequencing across eight internode segments. Previously described topologically associated chromatin domain borders aligned with transcription start sites of constitutive expressed genes. CpG and CHG differential methylation was enriched in young segments, consistent with maintenance methylation; whereas CHH methylation showed similar differential abundance in young and old segments. CHH and CHG methylation in upstream regions, CpG methylation within gene bodies, and any methylation in 5' and 3' untranslated regions were permissive of moderate to high gene expression. Very low expression was associated with CpG methylation upstream, CHG and CHH methylation within gene bodies, and CpG and CHG methylation downstream. A nonrandom subset of genes, including cell wall-related glycoside hydrolases, lignin and tricin biosynthesis enzymes, and WD40 proteins, showed methylation-expression correlations, with expression changes enriched in triple-marked elements. These results suggest that internode phenotypes of DNA methylation machinery mutants relate to alteration of specific target genes, opening approaches for grass culm improvement for lodging resistance and biomass production.

O_LIThe control of amount of nuclear proteins is fundamental in regulating plant gene expression, but the mechanisms of quantitative dynamics of the nuclear proteome are largely unstudied during adaptive responses to pathogens. C_LIO_LIHighly specific labeling, enrichment and measurement of the nuclear proteome was performed using TurboID LC-MS of Arabidopsis thaliana leaves treated with the pathogen-associated molecular pattern (PAMP), flg22, and/or cycloheximide. The chosen experimental approach allowed discrimination of the effects of translation, nuclear protein import, trafficking of preexisting proteins, derepression, and nuclear protein turn-over upon elicitation of basal immunity. C_LIO_LIThe highly specific, deep coverage of proteins in the nucleus makes this study a resource for anyone interested in plant nuclear proteome dynamics and defense. C_LIO_LIAround 2,000 nuclear proteins were repeatedly quantified, including more than 300 transcription factors or other proteins related to transcription. Several proteins with documented activity in endosomes were newly synthesized and imported into nuclei upon PAMP challenge, suggesting alternative nuclear functions in PAMP-triggered immunity (PTI). Circadian clock components, including the transcription factor, CIRCADIAN CLOCK ASSOCIATED 1 (CCA1)-HIKING EXPEDITION (CHE), were depleted upon PAMP challenge, suggesting a safeguard against untimely induction of systemic acquired resistance (SAR). C_LIO_LIBased on proteomic patterns, proteins moonlighting in the nucleus as well as trafficking and turn-over regulation of the proteome are common elements during plant immunity. C_LI

Most proteins targeted to the organelles of endosymbiotic origin are encoded in the nuclear genome, placing them under the regulatory dominance of the nucleus. For photosynthetic eukaryotes, nuclear-encoded chloroplast proteins arise via two routes: First, genes of cyanobacterial origin were relocated to the nucleus through endosymbiotic gene transfer (EGT). Second, proteins of eukaryotic origin emerged to support chloroplast function and structure. These proteins are reimported into the chloroplast via an import machinery. Reversing the transfer of such genes from the nucleus to the chloroplast genome may offer insights into chloroplast regulation and evolution. In this study, we established a highly efficient and accessible electroporation protocol for chloroplast transformation in the green alga Chlamydomonas reinhardtii, and used it to reverse-transfer two nuclear-encoded genes encoding proteins arising via the two routes described above: the cyanobacteria-derived chloroplast division protein FtsZ1 and the Rubisco-linker EPYC1 of eukaryotic origin. Regardless of origin, both chloroplast-encoded FtsZ1 and EPYC1 showed proper localization and functionality comparable to their nuclear-encoded counterparts. Together, our study provides a robust protocol for chloroplast transformation, a platform for investigating the evolutionary drivers of EGT, and a foundation for advancing chloroplast bioengineering. SIGNIFICANCE STATEMENTO_LIEndosymbiotic gene transfer has resulted in the mass migration of genes from the chloroplast genome to the nuclear genome. Reversing the gene transfer could reveal the evolutionary significance of genome partitioning. C_LIO_LIUsing the green alga Chlamydomonas reinhardtii, this study developed an efficient, electroporation-based protocol for chloroplast transformation. Relocating the genes encoding two chloroplast-targeted proteins, FTSZ1 and EPYC1, to the chloroplast genome showed that the proteins maintained normal localization and function. C_LIO_LIThe established transformation protocol facilitates systematic testing of reverse gene transfer to elucidate the potential evolutionary advantages of genome partitioning and opens new avenues for chloroplast bioengineering. C_LI

Copper (Cu) is an essential micronutrient whose bioavailability is strongly affected by soil physicochemical properties. During evolution, plants have developed mechanisms to flexibly adjust their metabolism to Cu status. Superoxide dismutases (SODs), including Cu/ZnSOD1 (CSD1) and FeSOD1 (FSD1), are key antioxidant enzymes regulated in Cu dependent manner in Arabidopsis thaliana. Examination of CSD1 cellular distribution and activity revealed that CSD1 is a nuclear and cytosolic SOD whose abundance and activity respond to Cu availability inversely to FSD1. Combined microscopic and biochemical analyses of Cu-dependent dynamics revealed that, unlike FSD1, CSD1 localization in guard cells (GCs) remains independent of Cu availability. CSD1 escapes miR398-mediated regulation in GCs through a cell type-specific splice variant (CSD1.2) that carries an altered miR398-binding site. In silico analyses indicate that this mechanism is also present in crop species. Functionally, the csd1 mutant showed reduced sensitivity to abscisic acid (ABA)-induced stomatal closure, a phenotype rescued by reintroducing CSD1. Biochemical and reactive oxygen species (ROS) level analyses indicate that CSD1.2 most likely acts independently of its canonical enzymatic activity in GCs and functions upstream of the ROS burst in the ABA signaling pathway. Together, we present a novel, cell-type-specific mechanism that safeguards ABA-driven stomatal closure under fluctuating Cu supply.

Protoplast isolation is widely used for plant functional genomics and single-cell analyses, but its impact on transcriptional and cell state dynamics remains incompletely understood. Here, we generated time-course RNA-seq data from leaf protoplasts of Arabidopsis, maize, and poplar, sampling at multiple time points following isolation, to systematically characterize global transcriptional dynamics across species. We identified two major drivers of transcriptional variation: a persistent protoplast isolation effect and a progressive time-dependent transcriptional program, which can be divided into early, middle, and late stages corresponding to an immediate stress response, metabolic and chromatin regulation dynamics, and sustained metabolic and proteostasis regulation, together with species-specific differences across stages. We observed a rapid loss of cell-type-specific transcriptional signatures within 6 hours in Arabidopsis and maize, whereas poplar showed a slower decline. Single-nucleus RNA-seq at 6 hours in maize confirmed attenuation of cell-type-specific transcriptional structure. Furthermore, leveraging this time-course dataset enables the identification of aberrant cell states in single-cell RNA-seq data, exemplified by clusters showing elevated activity of protoplast isolation-associated, middle-, and late-stage transcriptional programs characteristic of stress-like states. Together, our results provide a cross-species framework for dissecting protoplast-induced transcriptional and cell state dynamics and facilitate the systematic identification of stress-associated cell states in single-cell transcriptomic data.

O_LIChloroplast division is governed by a conserved protein machinery, yet empirical characterization of these regulators remains limited in rice, a primary target for C4 engineering. Increased chloroplast occupancy in bundle sheath cells is a hallmark of the C4 pathway and so manipulating division is a potential strategy to achieve this goal. C_LIO_LIThrough developmental transcript profiling and image analysis, we identified a discrete window of active chloroplast proliferation in rice leaves, coinciding with peak expression of conserved plastid division genes. Functional characterization via overexpression revealed regulatory behaviours distinct from those in Arabidopsis thaliana. Overexpression of OsFtsZ1&2 resulted in fewer, enlarged chloroplasts per bundle sheath cell, whereas OsMCD1&OsMinE restricted plastid expansion without altering division rates. Conversely, overexpressing OsPDV1&2 or OsARC6&OsDRP5B increased plastid size without affecting total count. When OsPDV1&2 were co-expressed with transcriptional regulator ZmG2, we observed modest increases in chloroplast size alongside reduced stomatal aperture, increased stomatal density, and higher intrinsic water-use efficiency. C_LIO_LIThe results define the temporal landscape of plastid biogenesis in rice and demonstrate divergence across lineages. Our findings suggest that manipulating the division apparatus is insufficient to drive C4-like chloroplast biogenesis in the rice bundle sheath, highlighting the complexity of plastid-host cell coordination in cereals. C_LI

Reconstitution of functional nitrogenase in plants requires the coordinated expression of the [Fe-S] cluster assembly proteins NifU and NifS. However, the extent to which these proteins interact with endogenous Fe-S metabolism and affect plant physiology remains unclear. Here, we compared NifU and NifS homologs from diverse diazotrophs to identify variants compatible with the plant chloroplast environment. Selected variants of Azotobacter vinelandii, Fischerella thermalis, and Marinobacter lutimaris were characterized by transient expression in Nicotiana benthamiana and stable transformation in rice. Plant-produced NifU was largely devoid of [Fe-S] clusters when isolated but retained strong capacity for in vitro [Fe-S] cluster reconstitution and apo-NifH activation in a Ft > Av >Ml gradient, indicating correct folding and function but limited cluster loading or stability in vivo. NifU and NifS expression in transgenic rice resulted in variant-dependent proteome and phenotype effects, with A. vinelandii-expressing lines exhibiting severe defects, F. thermalis lines showing intermediate phenotype, and M. lutimaris lines being indistinguishable from wild type. These results reveal a trade-off between the biochemical activity of NifU and NifS and their compatibility with host metabolism, which must be considered for successful nitrogenase engineering in plants. HighlightNifU/NifS homolog selection determines trade-offs between [Fe-S] cluster assembly activity and plant compatibility, identifying variants that minimize physiological disruption while supporting nitrogenase cofactor assembly in chloroplasts.

The sequence of leaf physiological thresholds underlying plant responses to water deficit is thought to be functionally coordinated; yet, to what extent this coordination is maintained across genotypes and environments remains poorly documented at the intraspecific level. We characterized the sequence of stomatal closure, turgor loss and xylem embolism in the leaves of two genotypes of the riparian species Populus nigra (DRA-038 vs. PG-31) subjected to control, additional nitrogen or additional potassium treatments. Under control conditions, embolism measurements using the optical vulnerability method showed that DRA-038 was more vulnerable than PG-31, in agreement with measurements performed on stems with the reference Cavitron method. Stomatal closure consistently preceded xylem embolism, while bulk leaf turgor loss was typically observed once xylem embolism had already reached 50%. Hydraulic thresholds responded to treatments in a genotype-dependent manner, the intrinsically more vulnerable genotype DRA-038 being typically more plastic. However, despite variations across genotypes and treatments, the trait sequence remained tightly coordinated such that stomatal safety margins (SSMs) remained virtually null. These findings support a strong mechanistic integration of leaf hydraulic thresholds in poplar across genetic units and varying environments, questioning whether to favour intrinsic tolerance or plastic capacities in breeding future drought-tolerant genotypes.

Flowering time strongly influences crop adaptation, plant architecture, generation turnover, and breeding efficiency, but the functional organization of florigen genes remains poorly resolved in many polyploid orphan crops. Quinoa (Chenopodium quinoa) is a climate-resilient allotetraploid crop with extensive variation in flowering behavior, and genome analyses have identified multiple FLOWERING LOCUS T (FT)-like homologs. However, genome sequence and expression information alone cannot determine which homologs provide effective florigenic output in planta. Here, we combined apple latent spherical virus-mediated overexpression (VOX) and virus-induced gene silencing (VIGS) in quinoa with heterologous expression in Arabidopsis thaliana, domain-swapping analyses, and cross-germplasm validation to functionally dissect quinoa FT activity. Although several CqFT homologs were transcriptionally induced during the floral transition, their functional outputs were markedly unequal. CqFT1A and CqFT1B-1 acted as the major florigenic activators: overexpression of either gene induced rapid and synchronized flowering, whereas CqFT1-VIGS delayed flowering. In contrast, CqFT2A and CqFT2B retained only weak flowering-promoting activity, whereas CqFT1B-2 showed no detectable promotive effect under the conditions tested, revealing a clear functional hierarchy among transcriptionally induced CqFT homologs. Domain-swapping analyses showed that C-terminal variation contributes to, but does not fully explain, functional divergence among CqFT homologs. In late-flowering highland lines, elevated FT input accelerated flowering, induced coordinated floral transition, and shortened the time to viable seed production. These findings identify CqFT1A and CqFT1B-1 as the major florigenic activators in quinoa and establish a functional genomics framework for resolving and modulating flowering-time control in polyploid orphan crops.

Adaptation of plant photosynthesis to dynamic light conditions experienced in natural environments is achieved through specific protective mechanisms. Energy-dependent non-photochemical quenching (qE), regulated by Photosystem II Subunit S (PsbS), is a key process facilitating acclimation to fluctuating light in C3 plants, which operate conventional photosynthesis. C4 plants, which include some of the world's most productive and agriculturally important crops, have evolved a distinct high-efficiency photosynthetic pathway. Little is known about the role of specific processes, like qE, in acclimation of C4 plants to dynamic light environments. We generated gene-edited lines of the model C4 grass Setaria viridis lacking PsbS, which were found to be deficient in qE. This deficiency resulted in a modest increase in PSII photoinhibition and a CO2 assimilation penalty under light stress in short-term experiments, but photosynthesis and growth under fluctuating light were unaffected. Instead, keeping Photosystem I oxidised through photosynthetic control, negative feedback regulation of the Cytochrome b6f complex, was critical. Therefore, unlike in C3 plants, qE does not provide a significant adaptive advantage to C4 plants under dynamic light conditions. These findings provide important insights into the biology of C4 plants and help prioritise future strategies for improving the productivity and resilience of C4 crops.

The continuous transfer of organellar DNA into the nuclear genome during eukaryotic evolution has resulted in the widespread occurrence of nuclear plastid DNA insertions (NUPTs) and nuclear mitochondrial DNA insertions (NUMTs). However, their functional significance in nuclear gene expression and genome evolution remains largely unresolved. In this study, we employed Oxford Nanopore Direct RNA Sequencing (DRS) to investigate the transcription of NUPTs and NUMTs in the Populus nuclear genome and compared their transcriptional characteristics with their genome-wide insertion patterns. Our analyses revealed that the majority of transcribed NUPTs and NUMTs are enriched within introns and are co-transcribed with their host or adjacent genes in polycistronic-like transcriptional units. In addition, NUPTs and NUMTs frequently generate intronless transcripts, features reminiscent of their prokaryotic ancestry. We further identified a putatively functional NUPT-derived psbH gene that is unique to P. trichocarpa, providing new insights into the evolution of nuclear-encoded organelle-targeted genes. In addition, we identified transcribed NUPT and NUMT insertion polymorphisms among alleles, suggesting that organellar DNA insertions contribute to allelic variation and may participate in environmental adaptation. Collectively, our findings reveal previously unrecognized roles of NUPT and NUMT transcription in gene regulation, allelic variation, genome evolution, and the emergence of novel genes.

Intraspecific hybridization in allopolyploid plant genomes has the potential to induce non-additive changes in gene expression and DNA cytosine methylation, partly through interactions among divergent parental subgenomes. However, the extent to which intraspecific hybridization reshapes gene expression, coordinates homoeolog regulation, and remodels methylation in higher-order polyploids remains poorly quantified. To address this, we sequenced seedling leaf transcriptomes and methylomes from two parental cultivars of hexaploid bread wheat (Triticum aestivum L.) and their hybrids. More than 40% of genes were differentially expressed between hybrids and parents, although many were not differentially expressed between the parents themselves, consistent with complex trans-regulatory effects in the hybrid genome. This effect was more pronounced for homoeologs whose relative expression differed between the parents. These expression shifts often occurred simultaneously across all three homoeologs within triads, reducing homoeolog expression bias (HEB) in the hybrids. CG methylation levels were similar between the parents and hybrids in regions of low genetic divergence and in transposable element (TE)-rich regions, whereas CG sites in gene-rich regions showed more additive inheritance (hybrids intermediate between parents), particularly when parental haplotypes were themselves divergent. TE and gene body methylation (gbM) was strongly conserved in parents and hybrids. gbM was associated with more balanced homoeolog expression and fewer non-additive expression changes. CHH methylation showed overdominance, whereas non-conserved CHG methylation was enriched in TE-rich regions, suggesting that non-CG remodeling may reflect parental differences in TE and small-RNA content. Our results show that intraspecific hybridization within a hexaploid species can generate non-additive changes in gene expression and DNA methylation in seedling leaf tissue, while the presence of homoeologous genes, parental HEB, parental genetic and methylation divergence, and genomic location have varying levels of influence on expression or methylation remodeling.